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Abstract 

Relativistic electron beam propagation in plasma is fraught with several micro instabilities like 
two stream, filamentation etc., in plasma. This results in severe limitation of the electron transport 
through a plasma medium. Recently, however, there has been an experimental demonstration of 
improved transport of Mega Ampere of electron currents (generated by the interaction of intense 
laser with solid target) in a carbon nanotube structured solid target [Phys. Rev Letts. 108, 235005 
(2012)]. This then suggests that the inhomogeneous plasma (created by the ionization of carbon 
nano tube structured target) helps in containing the growth of the beam plasma instabilities. This 
manuscript addresses this issue with the help of a detailed analytical study and simulations with 
the help of 2-D Particle - In - Cell code. The study conclusively demonstrates that the growth 
rate of the dominant instability in the 2-D geometry decreases when the plasma density is chosen 
to be inhomogeneous, provided the scale length l/k s of the inhomogeneous plasma is less than the 
typical plasma skin depth [c/uj p o) scale. At such small scale lengths channelization of currents are 
also observed in simulation. 
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I. INTRODUCTION 


The relativistic beam plasma system occurs in the context of many frontier research 
areas in both laboratory and space plasmas, such as the fast ignition scheme of inertial 


confinement fusion 


EL 


compact particle accelerators [) , solar flares physics |3j, cosmic 


fl, 


magnetic field generation 4] and magnetic field reconnection [5j. Therefore, the study of 
instabilities associated with the relativistic beam-plasma systems is of central importance in 
understanding issues pertaining to a better transport of beam electrons in the plasma and 
provide clues to improve upon it. The manuscript explores this in detail. 

When a short, intense laser pulse (I> 10 19 W/cm 2 ) irradiates a solid target, a relativistic 
electron beam is generated through the process of wave breaking j3, H, 1 7]. This relativistic 
electron beam can propagate inside the plasma target even when the current associated with 
it exceeds the Alfven current limit / = (mc 3 /e )7 = 1776 kA, here m is the electron mass, e 
is the electronic charge, c is speed of light and 7 ^ is the Lorentz factor of the beam. This 
is because simultaneously a return current is induced by the electrons of the background 
plasma which compensates the forward beam electron current. This combination of forward 
and return current is, however, highly unstable to two stream and filamentation instabilities. 
For the two stream (longitudinal mode) js] instability, the perturbations propagate parallel 
to the direction of the beam flow. When the perturbations propagate normal to the beam 
jropagation direction they are termed as the filamentation instability (transverse mode) 


10]. This is also often known as the Weibel instability fll ]. 


In the general case the perturbation would propagate in a direction oblique to the beam 
and would have both filamentation and two stream characteristics. The instabilities associ¬ 


ated with the obliquely propagating perturbations have been investigated in the co 


limit 


12] as well as by kinetic approaches 13]. The particle-in-ccll (PIC) simulations 


d fluid 

dm 


have been carried out to validate the results. It has been shown that when 7 & > 1, relativis¬ 
tic effects come in to play and the filamentation instability dominates over the two-stream 
instability 16]. A full hierarchical map showing the parameter regions (ratio of the beam 
density to the plasma density and 7 ^) where these modes dominate has been provided in 
|l7 1. The maximally growing mode amongst the spectrum of unstable modes determines 
the beam evolution and sets the stage for the subsequent non-linear evolution. In fact the 
magnetic field generation is also closely linked with these modes and occurs due to the cur- 
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rent imbalance created by the most dominant modes of the instability. As the amplitude 
of the perturbations increase the nonlinear coupling defines the spectral cascade to other 
scales making the magnetic power spectrum turbulent in many scenarios. This in turn is 
detrimental to beam transport. 

There are many attempts at avoiding and/or suppressing such instabilities with are re¬ 
sponsible for current separation. For instance, theoretical and experimental work on col¬ 
limated propagation of electron beam assisted by resistivity gradient have been reported 


in Ref. 


18 


19]. A recent experiment 20] has demonstrated transport of Mega Amperes 


of electron currents through nano structured target over distances of the order of millime¬ 
ter. The transportation of electron beam over such distances is possible only if the beam 
plasma associated instabilities are somehow suppressed and the separation between forward 
and return currents does not take place as rapidly as predicted by the growth rate of these 
instabilities. In a simplified 1-D treatment, the suppression of Weibel/hlamentation insta¬ 
bility by the equilibrium density inhomogeneity (due to the nanostructures) has been shown 
in some recent studies ]21, 22], However, for the parameters of experimental system 20], 


typically the oblique mode would be the fastest growing mode 


17]. This necessitates a 2-D 


treatment wherein the role of plasma inhomogeneity on the growth of oblique mode should 
be investigated. Such a study is the main objective of the work presented in this manuscript. 

We have carried out the linear analysis of the beam plasma instability in an inhomo¬ 
geneous plasma medium by using the coupled Maxwell and two-fluid electron equations in 
the relativistic regime. The inhomogeneity in the beam-plasma density has been chosen 
to have a sinusoidal profile orthogonal to the direction of the beam propagation. This in¬ 
homogeneous density mimics the scenario which would arise when a high power laser of 
intensity (I> 10 19 W /cm 2 ) incident on a structured target ionizes it. The propagation of 
the accelerated relativistic electron beam (REB) occurs at a very fast time scale at which 
the heavier ions are unable to respond and continue to maintain the inhomogeneous profile 
of their density. Such a inhomogeneous background ion density constrains the flow of the 
beam and background plasma electrons. The linear analysis of the coupled fluid Maxwell 
set of equation shows that when the plasma density ripples have a scale length less than the 
skin depth, the dominant oblique mode gets suppressed. We have also carried out a 2D PIC 
simulations which corroborates with the linear results. With the help of PIC simulations 
one has also been able to explore the nonlinear regime of the system. It is observed that 
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magnetic field energy is comparatively lower for the inhomogeneous cases in the nonlinear 
regime. This further suggests that the current separation is weak when the plasma is in¬ 
homogeneous, conducive to the uninhibited propagation of the REBs. The manuscript has 
been organized as follows. The theoretical model (coupled Maxwell and two electron fluid) 
and governing equation for linear analysis for the beam -plasma system has been presented 
in section II. The growth rate of the unstable mode has been analyzed in section III. In 
section IV, the observations from PIC simulation have been presented. A comparison be¬ 
tween homogeneous and inhomogeneous density plasma is also provided which shows that 
if the scale length of the inhomogeneity is sharper than the electron skin depth the growth 
of beam plasma instability gets suppressed. In section V, we present our conclusions. 


II. THEORETICAL MODEL AND GOVERNING EQUATIONS 


We employ a two fluid description wherein the relativistic electron beam (REB) and the 
background electrons are treated as separate fluids. The evolution equations of these two 
electron species are coupled with the Maxwell’s equations which governs the held evolu¬ 
tion. The ions are assumed to be static and merely provide a neutralizing background. A 
collision-less, quasi-neutral and un-magnetized electron beam-plasma system is considered. 
The electron beam drift in the y direction inducing a return current in cold background 
plasma electrons (T 0p ). The initial configuration is chosen so as to have charge and current 
neutrality. The sketch of the configuration is shown in Fig. [0 At time t=0, a space de¬ 
pendent (transverse to the beam how direction) background ion density rioi is chosen. The 
charge neutrality condition is fulfilled by choosing the sum of the beam density nob and the 
background electron density rio p to match with background density of ions. The system can 
be described by following set of dimensionless governing equations: 


dn a 

dt 


+ V ■ {n a v a ) = 0 


dfa 

dt 

dB 

~dt 

dE 

~dt 


+ iT a ■ Vp„ = - (e + i? a x B j 
= —V x E 
= V x B - E a J a 


VP Q 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


4 





with v a = Pa/(l+p« 2 )5 where p* a is momentum vector, J a = —n Q ;iT a . The pressure P a 
is provided by the equation of state. In the above equations, velocity is normalized by 
the speed of light c, density by n 0 (the spatially average ion density), frequency by ojq = 
\JAnnie/m e , length by electron skin depth d e = c/oj^e and electric and magnetic held by 
Eq = B 0 = m e cuJo/e where m e is electron rest mass and e is electron charge. The subscript 
a is 'b' for beam and 'p' for plasma. 

In the equilibrium there is no electric and magnetic held, so there is complete charge as 
well as current neutralization. This is achieved by balancing the forward and return electron 
currents at each spatial location. For simplicity of treatment, the prohle of transverse 
temperature T a is chosen in such a way that gradient of pressure is zero in equilibrium. The 
temperature parallel to beam propagation direction has been chosen to be negligible. The 
background plasma is chosen to be cold ( Tq p = 0) in all our analytical as well simulation 
studies. Thus, for an inhomogeneous case considered by us in this work we have the following 
conditions for equilibrium: 

riQiix) = n 0 b(x) + n 0p (x ) (5) 

^ a n 0a (x)voa = 0 ( 6 ) 

In equilibrium beam pressure P, b& is chosen to be independent of x. This is achieved by 
choosing the beam temperature T 0b (x) to be satisfying the following condition 

P ob = T 0b (x)n 0 b(x) = constant = C. (7) 

T ob {x) = C/n ob (x ) (8) 

The suffix 0 indicates the equilibrium fields. We linearize equations (l)-(4) to obtain linear 
growth rate of the instability. The inhomogeneity being periodic along 'x' we choose all the 
perturbed quantities to have the form of f a = ’ ±2 ' ± 1-0 / QJ e* ^ kx +Jk d x + k vy-^t). jq ere j s 
the wavenumber associated with the periodic prohle of inhomogeneity. For small amplitude 
of inhomogeneity (e < 1), the terms corresponding to j = 0, ±1 are the ones which are only 
retained. All higher order terms are neglected. The coupled set of differential equations 
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obtained after linearizing Eqs.(l)-(4) are following : 


n a7o a V "ay + ~ k y) ^«7o a^lay + ~ fcj) ^ V l»y + ^ K aVlay)' = 0 (9) 


(“ 2 - k 2 y)^, , f^ 2 - k D 


a 


"7-^0 a^’lax 


a 


VJ ^ ' ‘ + O 2 - &y)^a70a ) V lc 


V^O oft lax 




a 


-a,y noaKioi + (*7 (w 2 - KplTru. — i {k y — v„ a io) SVriL) >i i;/ = 0 


( 10 ) 


Where f2 a =(o; — k y voa ) and 77 is ratio of specific heat. When the plasma density is homo¬ 
geneous i.e. e = 0 , these coupled equations reduce to the standard linear equations for the 
beam plasma system: 


fi a7 L V lay + (<*? - k l) ^«7o a Vlay + (w* ~ tf) ^ ^ Tl 0a v[ ay = 0(11) 


Vtr 


Vtr 


(a ; 2 - kl) 


q VToaVi ax “I - ky ) a^lax ^ ^ ^0 o/VIolx 

\L a 

+ ~ k l)V T 0 a - l {ky - V 0a uj) ^a7o a ^) V lay = 0(12) 


In the next section we study in detail the growth rate of the most unstable mode for the 
general case by using eq. m- 


III. ANALYTICAL STUDY 

The linear evolution of beam-plasma system is studied by choosing sinusoidal form of ion 
background density of the form no* = 7io(l + ecos{k s x )) where n 0 is constant normalized 
density, £ is inhomogeneity amplitude and k s = 27im/L x ( m is an integer and L x is the system 
dimension along x) is inhomogeneity wave vector. The equilibrium beam and plasma density 
profile have been taken as n 0 6 = /3&no(l + £fccos(/c s x)), n 0p = /3 p no(l + £ p cos(k s x)) to maintain 
quasi neutrality where (3 a = 1 is a fraction. For simplicity, we have chosen e b = e p = e 
here. 

To evaluate the linear growth rate of oblique mode driven instability, we expand eq. (HUl) 
and eq. (fill) and solve it by substituting for the Bloch wave function form of the perturbed 
fields in the periodic system. Retaining only the first order terms, as mentioned earlier we 
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obtain: 


fi «7oa ( A2 - k d ) V day + ^ ^ + ^ J2 


0.5£/3 a 

”q v ay 


+ * ffP^dVdo 


+i 22 iT" - 0 .5s/3 a kdV ax = 0 (13) 

' A 


^«70a (-^ 2 - &o) Wy + A 2 T^Vay + A 2 ^2 ~7r~ (Vday + Vaay ) + * TT^ ok ° 


_ /■ 

+* ^ ^O.Se^fco (v daa; + rw) = 0 (14) 

a 


fi «70a ( A2 - k T) V aay + A 2 ^ 


Pa 


q v aay 

a Lfy 


+ A 2 £ 117k,, 


“1“ ^ ^ ^ q fiot.k'aVaoLX 
\ L a 
a 

+i 22 2r®' 5£ P ak a Vax = 0 ( 15 ) 


0, 


ay 


\ \2 

-r]c a k 2 d - A 2 H« 7 0q v dax + Q a 22 PaVdax + 0.5 e———i]c a [k 2 d - k 0 (2k 0 - k s )]v ax 

' J ^ ^ OiPoL 

' a 

/A 2 \ A 2 

+0.5efl a y ^PgVgx + ( k y k d -———r]c ol - C a kl a ^ a k d I v day - 0.5£——r)c a kyk o v ay = 0 (16) 

\ ^^aPa. J aPa 

\ ^2 

V C ak Q A ^qj'YOq: j ^ckcc "l - E fiaVax + 0.5e———//C q, [fcg — k d (2k d + A: s )]ud a x + 

/ ^aPa 


a,A 


A 


A 2 

0.5£— — ?7Cq;[/Cq k a (2k a ^s)]^aax “1“ O.S-Sf^Q, ^ ^ /^g {Vdax ^aa.x) 

^ ^olPol 

OL 

/A 2 \ A 2 

+ ( kyk()— ——T]C a Ca'^aToa^'O j 0.5<E^— — T/Co,ky {kdVday “I - k a V aOL y ) 0 (17) 

\ ^^aPa J ^^czPcz 

A 2 


\ ^2 

~^]^CX.k a A ^g^Oc* j ^aOLX ^g E ficxVaoLX + 0.5e—— -iic a [k 2 a - k 0 (2k 0 + k s )]v ax 
i J ^^aPa 

/ A 2 \ ^2 

+0.5eQ a 22 P»Vax + ( k y k d ——r]c a - Cabalaa k a ) v aay - 0.5£——ric a k y koV ay = 0 (18) 

\ ^aPa J ^ olPol 


where A 2 = (cu 2 - fc 2 ), k d = (k x - k s ), k 0 = k x , k a = (k x + k s ) and C =(k y - v 0a uj). 

The map of the growth rate in the parameter space of k x and k y for the homogeneous 
system has been shown in Fig. [2] for n 0b /n 0p = 1/9 or n 0b /n 0e = 0.1, v ob = 0.9c, v 0p = —0.1c 
and To&x = 10 KeV. This map shows that for this particular set of parameters, oblique 
mode with(fc x , k y —( 1.5,1)) has a maximum growth rate. Also a narrow oblique strip of 
wavenumbers extending up to k x =oo is found to be unstable. 

When the plasma density in chosen to be inhomogeneous and the scale length of the 
inhomogeneity is longer than the skin depth, i.e. 27 t//c s > c/ujq the growth rate of oblique 
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mode driven instability is observed to increase with the inhomogeneity amplitude £, as shown 
in Fig. [3] and [4] The unstable modes continues to stay around the oblique patch in the 
k x vs. k y plane. However, the growth rate is spread over a wider k y domain compared to 
the homogeneous case. Another feature is the appearance of several maxima along k x in 
contrast to the single extrema for the homogeneous case. 

When the inhomogeneity scale length is sharper than the skin depth, i.e. for the case of 
27 r/k s < c/ujq, there is an overall reduction of the growth rates in the system in comparison 
to the homogeneous case. Furthermore, with increasing amplitude of density inhomogeneity 
also the growth rate reduces in this regime. This has been illustrated in the 2-D plot of the 
growth rate in Fig. ED and [t)J 

A detailed comparison of the growth rate Y gr of the maximally growing mode for various 
different values of k s and e has been provided in TABLE I. 

TABLE I 

The maximum growth rate of oblique mode driven instability evaluated analytically under 
the approximation of weak inhomogeneity amplitude as well as from PIC simulation. 


£ 

k s 

r gr ( max .) 

F p ic ( max .) 

0.0 

0.0 

0.1822 

0.1823 

0.1 

7T 

0.1857 

0.1827 

0.1 

2t r 

0.1811 

0.1818 

0.1 

37T 

0.1728 

0.1813 

0.2 

7r 

0.1884 

0.1830 

0.2 

37T 

0.1475 

0.1659 

0.4 

37r 

0.1352 

0.16 


It is clear from the results that the growth rate drops when the amplitude as well as the 
wave number of the plasma density inhomogeneity is increased. 


IV. PIC SIMULATION 

In this section, we present the results from the 2|D PIC simulations employed for the 
study of beam plasma instabilities for an inhomogeneous density plasma. The simulation 
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box is a Cartesian X x Y plane with periodic boundary condition for both the electromag¬ 
netic field and charged particles. The initial choice of configuration satisfies the equilibrium 
condition. The electron beam is chosen to propagate in the y direction with a relativistic 
velocity whose current is neutralized by the cold return shielding current from the back¬ 
ground plasma electron flowing in the opposite direction with a velocity v 0p . The respective 
densities are appropriately chosen for the current to be zero. Charge neutrality and a null 
current density are both ensured initially (t = 0). Thus, the system is field free and in 
equilibrium initially. 

The electron beam has also been chosen to have a finite temperature T ob . For the inho¬ 
mogeneous plasma the pressure contribution in such cases have been avoided by choosing 
T ob to be space dependent defined by eq. |8l The ions are kept at rest during the simulation. 
The uniform plasma density no e is taken as 10n c where n c = 1.1 x 10 21 cm -3 is the critical 
density for 1/jm wavelength of laser light. The area of the simulation box R is 30 x 15 
(c/cuo) 2 corresponding to 1500 x 750 cells where c/u 0 = d e = 5.0 x 10 -2 /un is the skin depth. 
The total number of particles is 45000000 each for both electrons and ions. To resolve the 
underlying physics at the scale which is smaller than the skin depth, we have chosen a grid 
size of 0.02d e . The time step is decided by the Courant condition. The time evolution 
of the box averaged field energy density for every component of E and B are recorded at 
each time step. The energy density is normalized with respect to (m e CLOo/ e) 2 and time is 
normalized with respect to the electron plasma frequency ujq. The results obtained from the 
simulation of the beam-plasma system for various amplitudes and scale lengths of inhomo¬ 
geneity have been compared with the results of the homogeneous beam-plasma case. The 
initial choices of the parameters for simulation are taken as n 0b /n 0p = 1/9 or n 0b /n 0e = 0.1, 
t’ob = 0.9c, vq p = —0.1c and T ob ± = 10 KeV. These parameters favor the growth of the 
oblique mode over the filamentation and two-stream modes. While considering the simula¬ 
tion of the inhomogeneous case, the inhomogeneity amplitudes e is varied from 0.1 to 0.2. 
The inhomogeneity wavenumber is chosen as k s =7T, 2 tt and 37r for each of these amplitudes. 

The perturbed field energies (magnetic as well as electric) are tracked in the simulation 
for the study of instabilities associated with the system. The rate of exponential growth 
of the perturbed energy provides for twice the growth rate of the fastest growing mode 
associated with the instability. The growth rates T p j C (max.), thus evaluated for different 
cases of simulations have been presented in one column of TABLE I. 
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The linear regime can be clearly identified from a log plot of total energy shown in (Fig. [TJ). 
The region of constant slope after an initial transient provides the growth rate of the fastest 
growing mode. The straight line alongside represents the slope obtained from the linear 
theory provided in section II. It can be observed that the agreement between simulation 
and the theoretical linear results are remarkably good. The comparison also shows that the 
reduction in the growth rate in the presence of inhomogeneity in plasma density for scales 
sharper than the skin depth. 

The separate evolution of electric field energy density and the magnetic field energy 
density for various simulation cases have been shown in Fig. [H] and Fig. [2] respectively. 
The energy associated with electric field is observed to be typically always higher than 
the energy in magnetic field. It is interesting to compare these perturbed energies in the 
nonlinear regime. While for the homogeneous case the energy continues to remain high, 
there is a perceptible drop in both electric and magnetic field the energies in the presence 
of inhomogeneity in the nonlinear regime. This has important significance as it suggests 
that the forward and reverse currents which got separated during the linear phase have a 
tendency to merge again for the inhomogeneous density when nonlinearity sets in the system. 

The snapshot of the spatial profile in the x — y plane of certain fields ( B s , rib and n p ) 
at a time cuo t =33, 53 and 90 have been shown in Figs. QHl Qj] and [[2] respectively. In 
these figures the first column corresponds to the homogeneous case, i.e. £ = 0.0 and second 
column to the inhomogeneous case with £ = 0.2 and k s = 37t (corresponding to density 
inhomogeneity scale lengths to be sharper than the skin depth). It is clear from these 
figures that the variations in these fields are in both 'x' as well as 'y' directions, confirming 
that the oblique mode continues to dominate the instability scene. However, it should also be 
noted that during the linear phase the magnetic field acquires structures which are extended 
and aligned along the beam flow direction of y. The structure size along the x direction is 
significantly shorter and is found to compare with the typical values of the skin depth. In the 
inhomogeneous case additional structures at the shorter inhomogeneity scale are observed 
to ride over those appearing at the skin depth scale. As an aside we wish to mention 
that fluid simulations based on Electron Magnetohydrodynamic model has also shown that 
plasma density inhomogeneity leads to the guiding of electron current structures 23]. These 
studies, therefore, suggest that specially tailored targets incorporating appropriate forms of 
plasma density inhomogencity can in fact be helpful for efficient transport of electron beam 
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through plasmas. 


V. SUMMARY 

In this work instability of the beam plasma system has been analyzed for the case when 
the background plasma is chosen to be inhomogeneous. The linear analytical evaluation of 
growth rate shows that when the inhomogeneity scale length is shorter than the typical skin 
depth scale of the plasma the growth rate of the instability is suppressed. A quantitative 
comparison of the growth rate of the maximally growing mode for the homogeneous and 
inhomogeneous cases were made, which clearly show the reduction in the growth rate in the 
presence of density inhomogeneity with sharper scales. A detailed PIC simulation in 2-D 
have also been performed which confirms this. The simulation shows reduced growth rate 
and also the channelizing of current during the linear phase. 

Our work has direct relevance to the recent experimental work on the observation of 
efficient transport of Mega Ampere of electron currents through aligned carbon nanotube 
arrays. The ionization of the carbon nanotubes by the prepulse of laser pulse would produce 
an inhomogeneous plasma density. The suppression of the beam plasma instabilities would 
then aid the process of efficient transport of electron current as observed in the experiment. 
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FIG. 2: Growth rate map of the 2D oblique instability for a homogeneous beam-plasma at 

transverse beam temperature Xb)_L=10 keV. 
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FIG. 3: Growth rate map of the 2D oblique instability for a inhomogeneous beam-plasma 
at transverse beam temperature Tb 0_L = 10 keV at k s = tt and £ = 0.1. 
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FIG. 4: Growth rate map of the 2D oblique mode driven instability for a inhomogeneous 
beam-plasma at transverse beam temperature 71 0 _l=10 keV at k s = tt and e = 0.2. 
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FIG. 5: Growth rate map of the 2D oblique mode driven instability for a inhomogeneous 
beam-plasma at transverse beam temperature Tf,o±=10 keV at k s = 37t and e = 0.1. 
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FIG. 6: Growth rate map of the 2D oblique mode driven instability for a inhomogeneous 
beam-plasma at transverse beam temperature Tboj_=10 keV at k s = 37t and £ = 0.2. 
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FIG. 7: Growth rate from PIC simulation by energy slope. 
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FIG. 8: The evolution of y-component of normalized electric field energies Ue v with time 

for homogeneous and inhomogeneity cases. 
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FIG. 9: The evolution of z-component of normalized magnetic field energies Ub z with time 

for homogeneous and inhomogeneity cases. 
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FIG. 10: Time evolution and spatial configuration of z-component of magnetic field for 
homogeneous (first column) and inhomogeneous (e=0.2 at k s =37T , last column) case at 

time co 0 t = 33, co 0 t = 53 and u 0 t = 90. 
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FIG. 11: Time evolution and spatial configuration of beam density for homogeneous (first 
column)and inhomogeneous (e:=0.2 at k s =37T, last column) case at time oj 0 t = 33, u> 0 t = 53 

and u 0 t = 90. 
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FIG. 12: Time evolution and spatial configuration of background electron plasma density 
for homogeneous (first column)and inhomogeneous (£=0.2 at k s = 37t, last column) case at 

time ui 0 t = 33, co 0 t = 53 and oj 0 t = 90. 
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